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In this article we review the quantum Hall physics of graphene based two-dimensional electron 
systems, with a special focus on recent experimental and theoretical developments. We explain why 
graphene and bilayer graphene can be viewed respectively as J = 1 and J — 2 chiral two-dimensional 
electron gases (C2DEGs), and why this property frames their quantum Hall physics. The current 
status of experimental and theoretical work on the role of electron-electron interactions is reviewed 
at length with an emphasis on unresolved issues in the field, including assessing the role of disorder 
in current experimental results. Special attention is given to the interesting low magnetic field limit 
and to the relationship between quantum Hall effects and the spontaneous anomalous Hall effects 
that might occur in bilayer graphene systems in the absence of a magnetic field. 



I. INTRODUCTION 

The quantum Hall effect (QHE) is one of the most 
remarkable phenomena in condensed matter physics. 
From a phenonienological point of view it is a transport 
anomaly that is reinforced by disorder-induced localiza- 
tion of charged quasiparticles. It is observed in two- 
dimensional electron systems whenever a charge gap, i.e. 
a jump in the chemical potential, occurs at a density that 
is magnetic field dependent. The gap-density's magnetic 
field dependence can be directly related to Hall trans- 
port properties,'^ and to topological Chern indices that 
characterize the dependence of single electron wave func- 
tions on Bloch wavevector in the special case of periodic 
systems, and many electron wave f unc tions on boundary 
condition angles in the general case. 1^ The C hern i ndices 
determine the values of the Hall conductivitjff^EI] which 
are quantized in units of the conductance quantum e'^/h 
in non-interacting electron systems with fractional values 
allowed in interacting systems. 

The simplest quantum Hall gaps follow directly from 
the single-particle physics of kinetic energy quantization 
and Landau level formation in two-dimensional systems, 
but interaction induced gaps are also common. Lan- 
dau quantization gaps can occur only at integer values 
of the Landau level filling factor v. Interaction induced 
gaps can occur at integer or fractional values of v. In 
a system with no periodic potential, and therefore no 
length scal^i^ other than that defined by the magnetic 
flux quantum, gaps always occur at fixed u and have 
Hall conductance ve^ /h. Some of the incompressible 
states associat ed wit h interaction induced gaps have ex- 
otic propertiesj '^ l ^ l ^" ! including quasiparticles with frac- 
tional charge and fractional, or even non-Abelian quan- 
tum statistics. The discovery of the quantum Hall effect 
in the 1980's was a watershed event in modern physics. 
The observation of perfectly quantized transport coeffi- 
cients in a disordered system containing millions of elec- 
trons was not only a great example of nature's tendency 
to surprise, but also pointed the community tJl^ a new 



type of order - topological order - which can occur in 
condensed matter. 

Quantum Hall effects have been studied for several 
decades in the two dimensional electron gases (2DEGs) 
found at interfaces between different semiconductors and 
at interfaces between a semiconductor and an insula- 
tor. These systems continue to be actively investigatecP^ 
and to reveal surprises. Semiconductor heteroj unction 
2DEGs are not two-dimensional in the atomic sense since 
their electronic charge is typically spread over 30 or more 
atomic layers, but in the electronic sense because their 
transverse coordinate degree of freedom is frozen out 
by finite- size q uantization gaps. The recent isolation of 
graphenej ^^ l ^^ l a two-dimensional honeycomb lattice of 
carbon atoms, has enabled exploratiorfi^M! of the quan- 
tum Hall effect, along with many other electronic prop- 
erties, in a 2DEG that has a thickness of one atomic 
layer and is qualitatively distinct. Because one, two, and 
three layer graphene systems stacked in different ways 
have quite distinct electronic properties, graphene re- 
search has given rise to an entire family of 2DEG systems 
with a very rich set of behaviors. This review will focus 
for the most part on new questions that have arisen in 
experimental and theoretical studies of the quantum Hall 
effect in one and two layer graphene systems. 

The article is organized as follows. In section II 
we briefly summarize some of the main properties of 
graphene-based 2DEGs. In section HI we discuss the low- 
energy effective theories of graphene and bilayer graphene 
starting from the nearest neighbor tight-binding model. 
In the process we point out that the low-energy effective 
theories of single layer and bilayer graphene are partic- 
ular cases of a broader class which we refer to as chiral 
2DEG models. Other members of this class of Hamil- 
tonians are approximately realized in thicker few layer 
graphene systems. In section IV we analyze the proper- 
ties of chiral 2DEGs in the presence of a magnetic field, 
with special focus on their unusual Landau quantization 
pattern. In section V we focus on interaction driven 
quantum Hall effects in graphene and bilayer graphene 
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at integer filling factors, focusing on the unusual v = 
QH plateau of both monolayer and bilayer graphene. In 
section VI we address the recently discovered fractional 
quantum Hall (FQH) states in graphene and offer some 
speculations on the physics that determines which frac- 
tional states are strongest. Finally we conclude in section 
VII with some speculations on interesting directions for 
this rapidly evolving field. 

There are already a number of excellent review^i^'^oHllI 
on different aspects of the QHE in graphene, and there 
is bound to be some overlap of this review with others in 
the literature. The focus of this contribution is on the the 
Landau quantization pattern of graphene and few layer 
graphene chiral 2DEG systems, and on the rich family 
of interaction dominated quantum Hall states to which 
it leads. 



II. GRAPHENE QUANTUM HALL EFFECT 
PRIMER 

Graphene can be viewed as a single layer of carbon 
atoms that has been isolated from a bulk graphite crys- 
tal. Its electronic structur^^Sl consists of sp^ bonding 
and anti-bonding bands which are located far from the 
Fermi energy and irrelevant for transport, and two tt 
bands - one for each of the two 7r-orbitals per unit cell in 
the material's bipartite honeycomb lattice. In a neutral 
graphene sheet the lower energy tt valence band is full 
and the higher energy tt conduction band is empty. The 
gap between conduction and valence bands vanishes at 
two points in the Brillouin-zone, the inequivalent honey- 
comb lattice Brillouin-zone corners. Only states located 
at nearby momenta are relevant for transport. As de- 
scribed more fully in the next section, the quantum me- 
chanics of these states can be accurately described over 
a remarkably wide energy window using a k • p enve- 
lope function Hamiltonian which has the form of a two- 
dimensional massless Dirac equation in which the role of 
spin is played by the two- valued which sublattice degree of 
freedom. (Because there are two inequivalent momenta 
at which the gap vanishes, electrons also carry a two- 
valued valley degree-of-freedom when described using the 
continuum k • p theory in addition to their pseudospin 
and real spin degrees-of-freedom.) A number of inter- 
esting properties of 7r-electrons in graphene follow from 
this ultra-relativistic quantum mechanics, including the 
absence of a gap, linear dispersion at low energies, and 
momentum space Berry phases. 

In graphene. Landau level formation leads to un- 
usual Hall plateaus on which the Hall conductivity is 
quantizecpSl at half-odd-i nteger multiples (per valley and 
spin!) of e^/h. The ear Ijff^Iill observation of this prop- 
erty in graphene 2DEGs was an important demonstration 
that the experimental systems behaved nearly ideally as 
expected. When spin and valley degrees of freedom are 
included and Zeeman coupling is neglected, Landau level 
formation in graphene leads to quantized Hall conductiv- 



ities with values axij = ve^ /h with v = ±4(n -|- 1/2) = 
±2, ±6, ±10, • • • . The strongest quantum Hall effects in 
graphene occur at these filling factors. Linear band dis- 
persion in the absence of a magnetic field leads to Landau 
levels with energies (measured from the band-crossing 
Dirac point) that are proportional to N (see Fig. [l]). 
For this reason the Landau levels close to the Dirac point 
are more widely spaced allowing the quantum Hall ef- 
fect at v — ±2 to persist to room temperatures at very 
strong magnetic fields.^ Because of the availability of 
graphene 2DEG systems, the quantum Hall effect is no 
longer solely a low-temperature physics topic! Similarly 
the quantum Hall effect dX v = ±2 has been observed at 
low temperatures at magnetic fields strengths well below 
1 Tesla. 

As sample qualities improved experiments began to 
reveal a large number of additional Hall plateaua^^ not 
expected from Landau quantization alone, the plateaus 
associated with these gaps could only be induced by 
electron-electron interactions. Historically these gaps 
had been theoretically anticipated from two different 
directions. The take off point for one group of re- 
searcherd25H29| j^^d been the tendency, already familiar 
from semiconductor 2DEG work, toward gaps at inter- 
mediate integer filling factors whenever a group of Lan- 
dau levels are degenerate either accidentally or because 
of a Hamiltonian symmetry. The tendency for interac- 
tions in this case to generate mass terms in quasiparticle 
Hamiltonians that open gaps, breaking symmetries if nec- 
essary, can be understood in terms of the degeneracy of 
the Landau levels and the role of exchange interactions.!^ 
A common example of this type of inter action induced 
gap is spontaneous spin-polarizatiorpi^ that is accom- 
panied by a large exchange gap between Landau levels. 
(The Zeeman splitting is often very weak compared to the 
typical electron-electron interaction scale so that these 
gaps are essentially spontaneous.) For this reason the in- 
compressible states associated with gaps that lift Landau 
level degeneracies are often referred to as quantum Hall 
ferromagnets, even when the gap i n ques tion is not a spin 
gap. Another group of researcher^SIH anticipated gaps 
at filling factor = in graphene sheets based on the 
idea, referred to as magnetic catalysis, that these spon- 
taneous mass gaps nearly occur in the absence of a mag- 
netic field and are then catalyzed by Dirac-point Landau 
level degeneracies. Recent experiments seem to favor the 
quantum Hall ferromagnetism scenaricP^'f^ since inter- 
action induced gaps have been observed at many integer 
filling factors, not only at i/ = 0, in both single-layer 
and bilayer graphene. Interaction induced gaps were 
also anticipated at fractional Landau level filling factors, 
the fractional quantum Hall effect (FQHE) , and have re- 
cently been revealed"^^ '^^ by transport measurements in 
suspended graphene samples and in graphene on hexag- 
onal Boron Nitrid^^ (h-BN) substrates. 

Like single-layer graphene (see Figj2]for lattice struc- 
ture), Bernal stacked bilayer graphene is a gapless semi- 
conductor. (The most common structure of bulk graphite 
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FIG. 1. Schematic illustration of the integer QH efTect found in (a) conventional 2DEGs, (b) monolayer graphene, (c) bilayer 
graphene and (d) ABC stacked trilayer graphene as a function of chemical potential ^ measured from the carrier neutrality 
point for fixed magnetic field. The Hall conductivity a^y is shown in red (solid) line and the longitudinal conductivity in blue. 
Plateaus in these C2DEG systems appear at 4(n + J /2)e^ /h integer multiples of /h, where n is an integer and J determines 
the 2DEG chirality (see text for details). For monolayer, bilayer, and ABC stacked trilayer graphene J is equal to the number of 
layers. The chemical potential ^ is measured in units of h^/{ml%) for conventional 2DEGs and in units of {y/2hvD /IbY 

for a J-chiral 2DEG (see text for details). Note that Landau level separation increases with energy for J = 3, decreases with 
energy for J = 1 and is independent of energy for J = 2. The shaded blue region corresponds to the disorder broadened Landau 
level density of states, which is inverted for holes carriers to facilitate clarity. The disorder broadened Landau level density 
of states is roughly the same form as the longitudinal conductivity except that its peaks are typically less sharp. Because 
graphene is a gapless atomically thin 2DEG system its carrier density can be varied over a large range surrounding the carrier 
neutrality point. 



consists of stacked Bernal bilayers.) As shown in Fig. [3] 
Bernal bilayers are characterized by a relative displace- 
ment between adjacent honeycomb lattices, which places 
one of the sublattices in the top layer directly above the 
centers of the honeycombs of the layer below, whereas the 
other sublattice lies directly above a lower layer sublat- 
tice. Interlayer hopping between the closely spaced sites 
drives the energies of orbitals concentrated on the latter 
sites away from the Dirac point. The low energy degrees 
of freedom are associated with the remaining sites, the 
ones without a near-neighbor in the other layer. The low 
energy effective theory of bilayer graphene consequently 
also has a two- valued low-energy which sublattice degree- 
of-freedom, one which applies however only below the 
interlayer hopping energy scale ^ 0.3eV. (The model of 
a bilayer as consisting of two coupled two-component 2D 
Dirac systems applies up to a much higher energy scale 
~ 3eV.) Since the two low-energy sublattices are local- 
ized in different layers, in the bilayer case the model's 
pseudospin can also be viewed equally well as a which 
layer degree-of-freedom that is readily coupled to an ex- 
perimental probe simply by applying an electric field per- 
pendicular to the bilayer surface. 

The low energy conduction and valence bands in bi- 
layer graphene have parabolic dispersion^S- and a Berry 
phase of 27r, twice the value of Dirac fermions. In a 
magnetic field, bilayer graphene's Landau level quanti- 
zation leads to an other distinct integer quantum Hall 
effect (IQHEpS'^ with quantized Hall conductivities 



a.:^y = ve^/h a,t ly = ±4(n + 1) = ±4, ±8, ±12 • • • . The 
eightfold degeneracy of the zero-energy = Landau 
level is due to the inclusion of both n = and n = 1 or- 
bital states among the zero energy solutions of the Dirac 
equation. This additional degeneracy, along with valley 
and approximate spin degeneracy, leads to a set of eight 
nearly degenerate Landau levels which are ripe for quan- 
tum Hall ferromagnetism, the N = octet. Although 
investigations of electron-electron interaction effects in 
bilayer graphene are still in their infancy, recent experi- 
ments in high magnetic field have clearly demonstrated 
additional plateaua^^'^ due to electron-electron interac- 
tions.EH 

Because its TV = states in different valleys are lo- 
calized in different layers, bilayer graphene is partially 
analogous to double layer quantum Hall systems in semi- 
conductors with two quantum wells. These systems have 
exhibited unusual transport anomalies due to excitonic 
superfluidity'^ at total filling factor z/^ = I. In the 
case of graphene bilayers excitonic condensation is ex- 
pected to be richer due to additional Landau level de- 
gree of freedom.'*'* Exciton condensation, or equivalently 
spontaneous interlayer phase coherence, is likely to occur 
when any two Landau levels that are localized in different 
layers are degenerate. This type of order is sometimes re- 
ferred to as quantum Hall superfluidity and is a special 
case of quantum Hall ferromagnetism. It should occujl^ 
in bilayer graphene at = 1,3 when the electric field 
between layers is very small. 
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Few layer graphene systems stacked in different ways 
give rise to distinct electronic properties at relevant en- 
ergy scales which introduces an entirely new class of 
2DEG systems. Recent theoretical work has indi cated 
that except for very weak fields multilayer graphen^^^I^ 
systems can be classified as a new class of 2DEG systems 
from here on referred to as chiral 2DEGs (C2DEGs). 
C2DEG models provide an accurate description of the 
low-energy properties of few layer graphene systems with 
a variety of different stacking arrangements consistent 
with symmetries. Due to the unusual Landau quan- 
tization of the single particle Hamiltonian, IQHEs in 
C2DEGs are expected to be remarkably different from 
that of semiconducting 2DEGs. Two examples of this, 
which have already been discussed, are single layer and 
bilayer graphene which exhibit unusual IQHEs described 
by Dirac continuum models. At low-energies these sys- 
tems represent the J = 1 and J = 2 instances of the 
C2DEGd^ model. Apart from these very recently several 
groups have reported experimental studies on ABA and 
ABC stacked trilayer graphene ABC trilayer graphene 
systems are approximately described by a single J = 3 
C2DEG model. 




FIG. 2. Honeycomb lattice of a single layer graphite flake with 
one sublattice in red and the other sublattice in blue. The 
lattice points form hexagonal plaquettes with the basis vector 
indicated by a and b. In the continuum limit the sublattice 
degree of freedom may be regarded as a pseudospin. The 
reciprocal lattice Brillouin zone is also a hexagon with the 
Dirac points marked by K and K'. When momentum k is 
measured from the Dirac point at the K or K' Brillouin zone 
corners, conduction and valence band eigenstates have their 
pseudospins oriented in the ±p direction. This property is 
often referred to as chirality. The angle ipp is the momentum- 
dependent phase difference between wavefunction amplitudes 
on the two sublattices. In bilayer graphene the pseudospin 
orientation angle is 2ifp. 



III. LOW-ENERGY THEORY OF SINGLE- 
MULTI-LAYER GRAPHENE 



AND 



Graphene is a two-dimensional material with carbon 
atoms arranged at the vertices of a honeycomb lattice. 
A simple 7r-orbital tight-binding model for its electronic 
structur^31] suggests correctly that graphene is a gapless 
semic onduc tor with linear dispersion about two nodal 
point^i^'^ located at the inequivalent corners of a hexag- 
onal Brillouin Zone (BZ) K and K'. In an electrically 
neutral graphene sheet the nodal points coincide with 
the Fermi level. When carriers are present due to dop- 
ing or gating, the nodal points lie at the center of nearly 
circular Fermi lines. 

In the case of Bernal stacked bilayer graphene, the low- 
energy effective theory also yields nodal points at the BZ 
corners, but the dispersion is quadratic.!^ In this section 
we derive these low-energy theories for graphene and bi- 
layer graphene, emphasizing that they are special exam- 
ples of C2DEG models.ESl 



A. Low-energy theory of graphene 

The honeycomb lattice is a textbook example of a non- 
Bravias lattice, i.e. a lattice with a basis of more than 
one atom per unit celhl^It is illustrated in Fig. [2] where 
A and B sublattice sites are represented by red and blue 
circles. In this section, starting from a 7r-band near- 
est neighbor tight-binding description of electrons on a 
planar honeycomb lattice, we show that the low-energy 
states of graphene are equivalent to those of massless 
Dirac fermions. We choose the two basis vectors of a hon- 



eycomb lattice to be a = ax and b = — a/2x -I- \/3/2ay, 
where = 0.142nm is the carbon-carbon distance. 

Within the nearest neighbor tight-binding approxima- 
tion the Hamiltonian in the atomic site basis can be writ- 
ten as 





-70/* (k) 



^7o/(k) 




(1) 



where 70 (~ 2.8eV) is the parameter for nearest neighbor 
hopping between sublattices A and B and /(k) is defined 
as 



/(k) + + e'^-^^+'''>. 



(2) 



The TT orbital site energies are identical by symmetry and 
chosen as the zero of energy. The conduction and va- 
lence band energies i?(k) = ±7o|/(k)| both vanish at 
the two Brillouin-zone (BZ) corner points K and K'. 
These points are usually called Dirac points for rea- 
sons that will become clear towards the end of this sec- 
tion. The Dirac point positions are given explicitly by 
K = -27r/(3a)x - 27r/(y3a)y and K' = -K. Expand- 
ing /(k) around the Dirac point K with k = K + p gives: 

/(P) = P • |^|k=K = ^7oafe - iPa)- (3) 

The effective Hamiltonian for states that are close to the 
nodal K point is therefore 



H = vd<t ■ p. 



(4) 
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where are Pauli matrices that act on the sublat- 
tice degree-of-freedom, pi = —ihdi is an envelope func- 
tion momentum operator. Here vjj = (\/3/2)7oa ~ 
1 X lO^m/s is graphene's Fermi velocity which is approx- 
imately 1/300 times the speed of light. The sublattice 
representation for valley K' can be chosen so that the 
general Hamiltonian is 'H(p) = £,V]j(t ■ p with ^ = ±1 for 
valleys K and K' respectively. To be consistent the pseu- 
dospinor for valley K is defined as -0k — ('Pak'^bk) 
whereas the pseudospinor for valley K' is defined as 

V'k' = ('^B K'' "^Ii K') ^^'^ 4'a{4'b) are envelop wavefunc- 
tions for the A and B sublattices measured from their 
respective valleys. So the sublattice labels should be re- 
versed when interchanging the valleys. The band energy 
dispersion is linear near both Dirac points K and K': 
e(p) = ±vd\p\- 

The low energy band model for graphene, described 
above, is equivalent to a model for two-dimensional mass- 
less Dirac fermions with a material specific speed of light. 
When the Fermi energy lies at the Dirac points the sys- 
tem is neutral. The Fermi energy lies above (below) 
the Dirac point for electron (hole) doped graphene. The 
eigenstates around valley K are: 



i±,p; 



1 



K 



= -jVp/2 



y/2 V ±e^'^p/2 



(5) 



where (pp — tan^^ (py/p^) and -t- labels the conduction 
(tt*) band case and — the valence (tt) band. Note that 
the eigenstates change sign if the phase (p is rotated by 
2tt. This is the tt Berry phase mentioned earlier. 

Chirality for pseudospin 1/2 objects is defined as the 
projection of the pseudo(spin) operator in the direction 
of momentunJ^. 



a- p 



(6) 



From this definition of the chirality operator it is easy to 
see that it commutes with the Hamiltonian Q , and that 
|±,p)k and |±,p)k' are also eigenstates of the chiral- 
ity operator. Electrons (holes) have positive (negative) 
chirality for valley K (^ = 1). Chirality is a good quan- 
tum number only as long as Hamiltonian Q is valid, 
and therefore holds only as an asymptotic property that 
is precise when close to the Dirac points. The eigenstates 
around K and K' are related by time reversal symmetry 
leading to opposite chirality in the two valleys. 



B. Low-energy theory of Bilayer graphene 

The stacking arrangement of the two graphene layers 
in Bernal bilayers is illustrated in Fig|3] Like single-layer 
graphene, Bernal stacked bilayer graphene has gaps that 
vanish at K and K' and a neutral system Fermi level 
that coincides with the nodal energies. In this section we 




FIG. 3. Schematic illustration of the structure of bilayer 
graphene with two shifted honeycomb lattices. The direct 
hopping process between red sites in the bottom and top 
layer is represented by the hopping parameter 71. At the 
Dirac points this terms creates bonding and antibonding band 
eigenstates that are concentrated on the Ai and A2 orbitals. 
The low energy sites Bi and B2 are indicated in blue. The 
gapless conduction and valence bands have quadratic disper- 
sions and are concentrated on these lattice sites. Direct hop- 
ping between the low-energy sites (the 73 process) leads to 
trigonal warping of the band dispersions. 



start from the single-layer graphene low-energy Hamil- 
tonian which has four component spinors corresponding 
to the four carbon sites per unit cell in this 2D crystal. 
It is then possible to derive an effective two-component 
model which applies over a more limited energy range. 

The key to the bilayer graphene electronic structure 
is the spatial structure described in Fig. [3j We label the 
sublattice sites as Ai{Bi) and ^2(52) in the bottom(top) 
layers respectively. Two of the four inequivalent atomic 
sites {Ai and A2) are stacked on top of each other, while 
the others (i?i and B2) do not have a near-neighbor in 
the other layer. The Ai — A2 interlayer hopping am- 
plitude, ^Ai-A2 = 71 ~ OAeV, leads to the formation 
of a dimer state of the Ai — A2 orbitals. This process 
drives the energies of electrons on these sites away from 
the Fermi level, leaving one low-energy site per unit cell 
in each layer. An additional interlayer hopping process 
from Bi to i?2, IB1-B2 =73 << 7i leads to a trigonal 
warping of the band structure^' at the lowest energies. 
We neglect this and other remote weak interlayer hop- 
ping processes for the moment and will argue later that 
they can be neglected in the large magnetic field limit. 

We follow McCann and Falkc^^ by starting from a 
continuum model bilayer band Hamiltonian which adds 
sublattice-dependent interlayer hopping to decoupled 
2D-Dirac bands: 



St 



"^4x4 — 





71 



( vdP 
vdP 

71 VdP 

\ VDP J 



(7) 
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for valley K where p = Px + ipy and vd is graphene's 
Fermi velocity. The above Hamiltonian operates on a 
four-component envelope function spinor wavefunction: 

= (V'Li,K'V']ii,K'V']i2,K' V'k,K) where V denotes the 
envelope function on the corresponding atomic site. The 
energy spectrum (p) takes the form 



e^(p)=7iCOs(^)±^/«|,|pP+72eos2(^), (8) 



with a — 1,2. As expected from our earlier discus- 
sion, two of the four bands have energies (p) « 7i 
and €2 (p) ~ ^.nd are well separated from the Fermi 
level in systems with moderate carrier densities. The 
other two bands are degenerate at the Fermi level and 
exhibit quadratic dispersion for small p. Expanding the 
low-energy bands ei{p) and e^(p) for |e| ^ 71 gives 



6(p)=± 



2m' 



(9) 



resembling the quadratic dispersion of a massive particle 
with an effective mass m = 7i/2w|) w 0.054TOe, coinci- 
dentally a value that is not very different from the effec- 
tive mass in GaAs (to* = 0.064TOe). The range of validity 
of the two-band bilayer graphene model can be estimated 
more precisely using Eq.Q which implies the condition 
|e| ^ 71/4. We can therefore expect the two-band bilayer 
model, the J — 2 C2DEG model as we explain below, to 
be accurate up to energy ~ lOOmeV. The range of va- 
lidity of the C2DEG model is therefore much more lim- 
ited in the bilayer than in the single-layer graphene case, 
for which it applies up to the eV energy scale. The low- 
energy chiral behavior has already been observed in quan- 
tum Hall experiments on bilayer graphene.l^ We leave 
discussion of the Landau level spectrum to the next sec- 
tion. 

The low-energy effective Hamiltonian for Bernal 
stacked bilayer graphene can be derived using degener- 
ate state perturbation theory. The dimer states formed 
by vertical interlayer hopping between Ai and A2 form 
a high-energy space selected by the projection opera- 
tor Vh — diag(0, 1, 1, 0). The low-energy subspace is 
spanned by sites Bi and B2 selected by the projection 
operator = diag(l, 0, 0, 1). We proceed by decompos- 
ing as: 



(10) 



where Hq = VhT^'Ph and V = "H — "Hq. The first-order 
contribution in perturbation theory vanishes leaving us 
with a second-order contribution to the low-energy effec- 
tive Hamiltonian: 



which for ii^ <C 71 gives 



1 

E-Uq 



-VhVVl, 



n 



(2) _ 
eff - 



71 







(pt)2 





(11) 



(12) 



^ (2) 

Here 'H Jt acts on the spinor containing only the low- 



energy sites = (V'bj K' k)- Identifying m 
ji/2vjj gives the same quadratic dispersion as before. 
The eigenstates for ( |l2| ) are obtained by substituting 
ipp 2(pp in Eq. 

In Bernal stacked bilayer graphene the eigenstates ex- 
hibit a Berry's phase of 2tt while the conduction and 
valence bands exhibit gapless parabolic dispersion. This 
2DEG system is qualitatively distinct from both single- 
layer graphene and normal semiconducting 2DEGs. In 
the next section we introduce effective Hamiltonians for 
C2DEGs, which include the low-energy effective Hamilto- 
nians of graphene and bilayer graphene as specific exam- 
ples. C2DEG models can be distinguished by a chirality 
index J. Graphene and Bernal stacked bilayer graphene 
are members of this class with chirality indices J = 1 



/,t 



and J — 2 respectively. In ( 12 1 above the subscript 2 can 



refer either to the order of perturbation theory or to the 
chirality index J — 2. 

Including the remote interlayer hopping 73 process and 
allowing for the possibility of an onsite energy difference 
between layers Ay gives additional contributions to the 



effective Hamiltonian:-'' Ti. 



38 1/(2) 



eff 



A, 



p 

pt 



2m7i 





n 



p^p 



(2) 
eff 



hyj + Ha where: 



(13) 



2m7: 



■pp 



where V3 = (\/3/2)73a. 73 leads to a distortion of the 
band structure that is invariant under 120° rotations, 
trigonal warping, which does not play an essential role 
in the high magnetic field limit. The external potential 
difference Ay between layers opens a gap in the spec- 
trum and modifies the shape of the conduction and va- 
lence bands. The modified dispersion has a Mexican hat 
shap^^ in which the band extrema occur along circular 
lines. 

The effect of interactions on undoped bilayer 
graphene's unusual band structure, which we have ne- 
glected to this point, is particularly interesting. Due 
to the absence of a Fermi surface, interaction effects at 
the charge neutrality point can lead to non-Fermi liq- 
uid behavior with ill defined quasiparticle excitations and 
a log squared behavior of the self-energy.l^ It has also 
been suggeste d tha t interactions can lead to the forma- 
tion of gappecP^*^ or gaplesJ^ broken symmetry ground 
states. Recent experiments on suspended dual gated bi- 
layer graphene samples do indicate anomalies near zero 
carrier densityin bilayer graphene in the absence of a 
magnetic fielcP^HSO] -^yj-^ich we discuss further below. 



C. Two Dimensional Chiral Fermions 

Graphene and Bernal stacked bilayer graphene are 
unique and distinct two-dimensional materials whose 
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low-energy electronic properties are well described by chi- 
ral fermion models.'*^ The single particle chiral Hamilto- 
nian with a chirality index J is 



e'ipy 

('ipy 



(14) 



where aj has the dimensionality of e/|p|'^, p = Px + iPy is 
momentum measured from nodal points in the Brillouin- 
zonc (BZ), and ^ = ±1 labels the two nodal points K 
and K' in the BZ. The pseudospin doublets in the re- 
spective valleys K(^ = -1-1) and K'(^ = —1) are defined 



as $ 



At J,t 



). From the 



analyses of the earlier sections we can see that graphene 
and Bernal stacked bilayer graphene's effective Hamilto- 
nians belong to this class with a chirality indices J — 1 
and J — 2 respectively. The higher J models apply ap- 
proximately to ABC stacked J-layer systems; although 
the relative importance of corrections to the model in- 
creases with J and the energy range over which it ap- 
plies approximately tends to decrease with J. For gen- 
eral N-layer graphene systems, the Hamiltonian can be 
approximately decomposed at low energies into a set of 
decoupled C2DEG models with J^i Ji = N and individ- 
ual Ji that depend on the detailed stacking sequence.'^ 
The energy dispersion of (14) is ej(p) — ztajp"^ with 
chiral band eigenstates: 



l±,P) 



1 / 

71 ±e^'^'^p/2 



(15) 



Just as before the pseudospin spinors in the different 
valleys have opposite chirality. 

These systems can be identified as different classes 
of zero-gap semiconductors. In each case the positive 
(negative) energies are identified with the conduction 
(valance) bands. The energy bands exhibit degenerate 
points K and K' in momentum space which are defined 
as the points where the conduction and valence bands 
meet. In neutral structures (i.e undoped with holes or 
electrons) the Fermi energy lies at the these degenerate 
points. 

Chiral Hamiltonians satisfy time reversal and inversion 
symmetries. The presence of two valleys is crucial for the 
time reversal invariance of chiral Hamiltonians. Time 
reversal symmetry for a chiral Hamiltonian refers to the 
property that (H (g) ax)'Hj{p){n ^ cTx) = Hj(— p) where 
n swaps ^ = -1-1 and ^ = — 1 in valley space. Chiral 
Hamiltonians also have the spatial inversion symmetry 
property (H ® I)-H,7(p)(H (g)I) = 'Hj{-p). 

The quasiparticles described hy Hj acquire a Berry 
phase of Jn upon adiabatic propagation along a closed 
orbit. This has unusual consequences for the single 
particle properties of chiral systems, the most notable 
of which is the unusual Quantum Hall effects that have 
been m easured in transport experiments in some of these 
systems.'^^'^^'^ Graphene and Bernal stacked bilayer 
graphene remain the primary focus of this review, 
although some important properties have systematic J 



dependences which we shall mention next. 



IV. THE INTEGER QUANTUM HALL EFFECT 
OF CHIRAL FERMIONS 

In this section we evaluate the spectra of chiral fermion 
models in the presence of a magnetic field, and discuss the 
unusual IQHEs which follow from these spectra. We also 
review transport experiments which verify the applicabil- 
ity of these models to single layer and bilayer graphene. 



A. Chiral Fermions in a Magnetic Field 

In the presence of a uniform magnetic field applied in 
the direction perpendicular to the plane of the C2DEG 
(B = Bz), the momentum operator in (14 1 is modi- 



fied to TT = p + {e/c)A where A — Ax + iAy, and A 
is the vector potential with B = V x A. We define the 
usual Landau level raising and lowering operators: 
and a with = {Ib/V^K)tt^ where Ib = (hc/eBy^"^ 



25.6/(-\/-S[Tesla])nm is the magnetic length. Eq. (14) 
for ^ = +1 then becomes : 







a 




(16) 



One immediate property of this Hamiltonian is the ap- 
pearance of zero-energy eigenstates of the form — 
(0, 0„) for n = 0, • • • , J — 1. Here the 0„ is the n*'* 
Landau level wavefunction for an ordinary 2DEG which 
is annihilated by a*"' for any k > n. Note that these zero- 
energy eigenstates are localized on the f (i) pseudospin in 
the K(K') valley, which for graphene corresponds to lo- 
calization of these Landau levels on sublattice A (B) and 
for bilayer graphene to localization on the top (bottom) 
layer. This property leads to a 4J-fold degenerate (in- 
cluding valley and spin) zero-energy Landau level. This 
degeneracy is, of course, on top of the normal Landau 
level one- state-per- flux- quantum degeneracy. 

The presence of wavefunctions with spatial structures 
that appear at different energies in the ordinary non- 
relativistic 2DEG model in the same degenerate manifold 
can create some terminological confusion. We will refer 
to the wavefunctions 4>„ as Landau level n orbitals and 
(as already anticipated) use upper case letter N to dis- 
tinguish levels with different Landau-quantized band en- 
ergies in the C2DEG model in a magnetic field. Thus the 
iV = Landau level eigenstates of a C2DEG model in- 
cludes states with orbital Landau level indices from n — 
to n = J — 1. 

The energies of higher N Landau levels are given by: 
eN,j = saj(^^^ ,/N{N-l)...{N-J+l), (17) 
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where s = ± distinguishes positive energy conduction 
band (s = +1) and negative energy valence band (s = 
— 1) Landau levels. The corresponding eigenstates are 

for valley K. The eigenstates of these Landau levels are 
a linear superposition of 2DEG Landau level orbitals n 
on one sublattice and n — J on the other sublattice. 

The energy spectrum of a chiral system in a magnetic 
field is remarkably different from that of a normal 
2DEG in which the Landau levels are equally spaced 
with the gap between adjacent levels equal to hiOc 
{ujc = eB/mc), and no zero energy eigenstates. For 
a chirality J system the level spacing energy scale 
varies with field as B^^"^. Spacings are larger near 
the Dirac point for J — \ and smaller near the Dirac 
point for J > 2. The unusual Landau level struc- 
ture of single-layer and bilayer graphene has been 
directly measured in scanning tunneling microscopy 
experiments o n grap hite subjected to a perpendicular 
magnetic field .1^^^ The Landau level spectrum for 
chiral fermions leads to a family of anomalous IQHEs 
which have now been verified experimentally in single 
layer, bilayer, and trilayer graphene. In next subsection 
we discuss the unusual IQHE of chiral fermions explicitly. 



B. Integer Quantum Hall Effect 

As explained in the introduction, the quantum Hall 
effect is a transport anomaly which signals a jump in 
the chemical potential, a charge gap^ at a magnetic field- 
dependent density. In translationally invariant contin- 
uum models the only length (and hence density) scale is 
the magnetic length, which measures the area that en- 
closes one quantum of magnetic flux. It follows that 
charge gaps in uniform 2DEG systems can occur only 
at fixed values of the filling factor z/, which is defined as 
the number of electrons per flux quantum. (In graphene 
systems the electron density is measured relative to the 
density of a neutral sheet.) A charge gap at filling fac- 
tor ly, when combined with incompressible state quasi- 
particle localization, leads to plateaus along which the 
longitudinal conductivity axx and the Hall conduc- 
tivity axy — >■ ve^/h. The use of the word plateau em- 
phasizes that because of charged excitation localization 
in the limit T — > these values apply over a finite range 
of magnetic field, of carrier density, or of any other ex- 
perimentally controllable parameter on which the 2DEG 
Hamiltonian depends. The relationship of the Hall con- 
ductivity to Chern numbers can be used to show that the 
quantum Hall effect can occur only at integer values of 
in the absence of electron-electron interactions. 

Gaps at rational fractional values of i/ (usually ones 
with odd denominator^^ have been studied for decades 
in ordinary 2DEGs and are now starting to be studied 




n (1 0" erir^ 

FIG. 4. (Adapted from RefUZlsee also Ref,^ ) The Hall 
conductivity a^y and the longitudinal resistivity pxx of mono- 
layer graphene as a function of carrier density at i? = 14 T 
and T = 4 K. This data demonstrates the half-integer QHE 
with chiral index J — 1 (graphene). The inset shows the Hall 
conductivity a^y in bilayer graphene with a chirality index 
J = 2. 



in graphene systems. They can be produced only by 
many-body electron-electron interaction effects. Incom- 
pressibility at a magnetic field dependent density always 
leads to gapless excitations localized at the edges of the 
system, the smoking gun signature of the QHE.— 

In C2DEGs, as in any other system, incompressibility 
at a magnetic field dependent density implies Hall quan- 
tization. The filling factors at which gaps are produced 
by Landau-level formation alone differ however. When 
the Zeeman energy and Ay are set to 0, gaps appear at 
the following set of integer filling factors (as indicated in 
Fig. |4]): 

iy^±4{N+^), Ar = 0,l,2,... . (19) 

Note in particular that there is no gap at zero carrier 
density v = 0. The above quantization condition fol- 
lows directly from the energy levels of chiral fermions in 
the presence of a perpendicular magnetic field discussed 
above. 

To understand the edge state implications of this un- 
usual quantization condition let us first review the IQHE 
in a normal 2DEG in which Landau level quantization 
leads to the energies e„ = (n -|- l/2)huJc- Landau levels 
are equally spaced with spacing equal to the cyclotron 
energy hiOc- In the presence of disorder Landau levels 
get broadened. When the Fermi energy lies in a gap be- 
tween Landau levels, the only states at the Fermi energy 
when disorder is weak are chiral edge states which can 
carry current without dissipation, causing the longitu- 
dinal resistance to fall to zero.^ When n Landau levels 
are filled, each Landau level gives rise to a gapless edge 
state branch. Because these n channels carry current in 
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FIG. 5. Schematic representation of the energy dispersion of 
the edge states of a graphene strip in a magnetic field within 
the continuum model. The electron-like (hole-like) band dis- 
persions are indicated in blue (red). The |A'^| > electron 
and hole-like bands are four-fold degenerate (spin and valley), 
whereas the A'' = bands are two-fold spin degenerate, xl 
and xu represent the respective left and right sample edges. 



Experiments now clearly show that at higher fields and 
in higher quality samples additi onal integer QH pla teaus 
appear for both single layei^^^ISll and bilayer graphen^^SEl 
that are not predicted in the Integer Quantum Hall se- 
quence (see Eq. 19). These new plateaus cannot be ex- 
plained by Landau quantization alone. They must origi- 
nate from the effect of Coulomb interactions or from ex- 
plicit symmetry breaking one-body terms in the Hamil- 
tonian - for example in the case of bilayer graphene from 
a large electric potential difference between the two lay- 
ers. Zeeman splitting is not on its own strong enough to 
produce additional quantum Hall effects in samples with 
typical disorder strengths, and must be enhanced by in- 
teractions. These additional gaps at integer filling factors 
are normally due to interaction-driven broken symme- 
tries and have been anticipated based on both Q uantum 
Hall FerromagnetisnPUMland Magnetic Catalysij S3 | 34 | 65 | 
scenarios, as discussed earlier. In the next section we dis- 
cuss these additional integer plateaus at greater length. 



parallel, axy = ne^/h. Each time the chemical potential 
crosses a Landau level there is an extra contribution from 
a new gapless edge state channel and the Hall conductiv- 
ity jumps by e'^/h. 

The unusual edge state structure and edge transport 
unique to chiral fermions can be understood from the 
Landau level spectrum and, in particular, from the pres- 
ence of zero-energy eigenstates.^'^ We can identify the 
Landau levels of a chirality J C2DEG as electron-like for 

> and hole- like for A'^ < as indicated in Fig. [5] for 
graphene {J ~ 1). The 4J-fold degenerate N = zero- 
energy state is neither electron-like nor hole-like. Given 
the spectrum of Landau levels for chiral fermions, the 
Hall conductance (Jxy must jump by AJe^/h when the 
chemical potential crosses the N — Landau level. The 
4 X J/2 term in Eq. (19) accounts for the presence of 



4 J zero-energy Landau levels which give rise to 4 J edge 
state branches. The unusual property of the N = Lan- 
dau level is that, as what one might guess on the basis 
of particle-hole symmetry, half of its edge branches ap- 
pear at energies above the bulk Landau level and half at 
energies below. The Hall conductivity on each plateau is 
therefore equal to the number of edge state branches in 
the corresponding gap. 

Soon after the isolation and identification of graphene 
this unusual IQHE was measured i n Hal l conductiv- 
ity experiments in both single layeJ^^^^ and bilayer 
graphen^^. Fig. |4] shows the experimentally mea- 
sured Hall conductivity for monolayer and Bernal bilayer 
graphencfSSl For single-layer graphene the gap between 
TV = and \N\ = 1 Landau levels is around 1300 K at 10 
T compared to around 100 K in an ordinary 2DEG. This 
is the reason the quantum Hall effect can be observed 
at room temper at urej^^S a great surprise given historical 
norms. 



V. INTERACTION INDUCED INTEGER 
QUANTUM HALL STATES 

In semiconductor 2DEG systems a wide variety of 
quantum Hall states have been discovered which are 
driven by interactions, i.e. incompressibilities develop at 
filling factors v at which the band Hamiltonian does not 
have a gap. From a fundamental physics point of view 
the most important of these are the ones that appear 
when a Landau level is fractionally filled, since they lead 
to fractionally charged excitations with fractional (and 
in some cases non-Abelian) statistics. Another impor- 
tant class of examples in semiconductor 2DEGs are the 
gaps which appear at odd integer filling factors. These 
are often found experimentally to be much larger than 
the Zeeman energy gap of the single-particle Hamilto- 
nian. In fact gaps do appear at odd integer values of 
v even when the Zeeman energy vanishes.— The ex- 
perimentally measured gaps are usually essentially pure 
exchange gaps due to electron-electron interactions and 
spontaneous spin polarization, i.e. due to what is of- 
ten referred to as quantum Hall ferromagnetism. The 
notion of quantum Hall ferromagnetism has been gen- 
eralized to other broken symmetry induced gaps in the 
quantum Hall regime.'^*' A particularly interesting aspect 
of quantum Hall ferromagnetism in semiconductors is the 
presence of skyr mion topological spin textures that carry 
electrical charge.^^i^ 

Band energy quantization for C2DEGs in the presence 

of magnetic field leads to a quantum Hall sequence with 

f = ±4{N + J/2) with V = ±2, ±6, ±10, • • • for single 

layer and i' — ±4, ±8, ±12, • • • for bilayer graphene. The 

four-fold degeneracy evident in the measured Hall con- 

duc tivity o f graphene is due to spin and valley degener- 
iVIls'lT " 



acy.t 



23] 



Bilayers exhibit an eight-fold degeneracy for 



iV = only which is due to the additional orbital de- 
generacy between orbitals with n = and n = 1 char- 
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acterPSSlwhen the spin Zeeman effect is inchided gaps 
appear at aU even integer filhng factors for single and 
bilayer graphene, except in the bilayer case for v ~ ±2. 
It was predicted theoretically that gaps would appear 
at all integer filling factors in samples with disorder en- 
ergy scales that are weaker than electron-electron interac- 
tion scales in both single layei^^SI and bilayei'^ graphene. 
This pr ediction has now been verified for sing le layer 
graphen^i^MIlll and also in bilayer graphen j^^ l ^^ l In this 
section we discuss the physics of interaction driven in- 
teger Quantum Hall states in graphene and in bilayer 
graphene. 

Single layer graphene gaps for v ^ ±4(A^ -f 1/2) can 
be due either to explicit symmetry breaking by one- 
body terms not included in our earlier analysis (that are 
usually enhanced by interaction dressing), or to purely 
spontaneous symmetry breaking by interactions. Inter- 
action driven additional integer quantum Hall states have 
been anticipated by Quantum Hall Ferromagnetisrr?^ '^ 
(QHF) theories which focus on interactions between 
sta tes in th e same Landau level, and by Magnetic Catal- 
yg^^^MB (MC) theories which focus on the v — case 
and assume that Landau level mixing plays an essential 
role. These theories can be further distinguished by their 
experimental predictions. While QHF theories predict 
that charge gaps and hence quantum Hall effects should 
be present at all integer filling factors, MC theories pre- 
dict that extra charge gaps should appear at = 0. The 
recent experimental observation of QH plateaus at all 
integer filling factorJ^ for graphene on h-BN seems to 
support QHF theories. 

Bilayer graphene also exhibits spontaneously broken 
symmetry states due to the presence of electron-electron 
interactions. In bilayer graphene QHF theories also pre- 
dict ordered states which exhibits incompressibilities at 
all integers plateaus. The N broken symmetry states 
are expected to be characterized by spontaneous spin 
and/or valley polarization. Because the = Lan- 
dau level includes an extra degree of freedom associated 
with the n = and n = 1 orbital degree of freedorrPsl 
the physics associated with its polarization, which is ex- 
pected at odd total filling factors, is more interesting. 
These orbitally broken symmetry states have unusual 
Goldstone collective modeJ^ associated with the dipoles 
induced by inter-orbital coherence. Apart from this the 
parabolic dispersion of bilayer graphene already makes 
it susceptible to interaction-induced symmetry breaking 
even at zero breaking magnetic field.'^^HlllExperiments on 
doubly gated suspended bilayer graphene samples indi- 
cate that incompressibilities exist for vanishing magnetic 
fields,'^ although the nature of the associated broken 
symmetry has not been unequivocally^^^ established as 
of this writing. 



A. Interaction induced integer QH states in single 
layer graphene 

1. Quantum Hall Ferromagnetism in Graphene 

Soon after the initial verification of the unusual IQHE 
of graphene, experiments^ at higher fields {B > 20T) 
revealed further plateaus 1^ = 0, ±1 and ly = ±4 that do 
not belong to the IQH sequence of graphene. These ad- 
ditional states suggest that the four-fold spin and valley 
degeneracy is lifted in a variety of filling factor-dependent 
broken symmetry states. Early tilted field experiments 
appeared to suggest that the states at i/ = 0, ±4 are 
spin polarized, implying that spin degeneracy is broken 
at these filling factors, whereas the v = ±.1 states depend 
only on the perpendicular component of the magnetic 
field, indicating that their gap is not spin related and 
corresponds instead to broken valley degeneracy.l^ How- 
ever, more recent experiments indicate that the v = Q 
gap is not spin-polarized and that the v — ±1 gap in- 
creases with in plane field. Temperature dependent ac- 
tivated measurements indicate that the gaps at = 0, ±4 
are larger than those expected from Zeeman splitting 
alone. If these gaps are spin-splittings, they are appar- 
ently strongly enhanced by many-body effects. Recent 
studies of graphene on h-BN substrates, which are much 
more ordered than Si02 substrates, appear to indicate 
that a number of the gaps in higher Landau levels in- 
volve both spin and valley degrees of freedom. 

Theories of quantum Hall ferromagnetism focus on in- 
teractions within a Landau level and the broken symme- 
try states they can induce. The discussion in the follow- 
ing paragraphs follows this line of thought and neglects 
Landau-level mixing effects. In these theories it is there- 
fore convenient to introduce a reduced filling factor for 
each Landau level N: 

m = ^<^ (20) 

where N^. is the number of electrons in the iV*'' Landau 
level and is the number of flux quanta enclosed in the 
system. In the continuum theory of graphene the spin 
and valley degrees of freedom are equivalent flavor labels. 
When Zeeman coupling and other weak perturbations are 
neglected, the Hamiltonian is SU(4) invariant. At integer 
fllling factors Quantum Hall ferromagnet ground states 
spontaneously break the SU(4) symmetry, ensuring good 
exchange energy for flavor indices within a Landau level 
and inducing a charge gap for excitations. 

The integer filling factor ground state of the SU(4) 
invariant Hamiltonian is known exactly. For filling factor 

i*o)= n ni.io)' (21) 

1 < o" < z> iv k 



is an exact eigenstate of any SU(4) invariant electron- 
electron interaction Hamiltonian. In ( 21 ) ^ is the elec- 
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tron creation operator for intra-Landau level orbital in- 
dex k and flavor (valley/spin) label a. |0) corresponds 
to the vacuum state in which Landau level N is com- 
pletely empty and all lower Landau levels are completely 
filled. The wavefunction |^'o) has lower energy because 
it puts many electrons into the same pseudospin state, 
minimizing their exchange energy to the maximum de- 
gree allowed by the Pauli exclusion principle at filling 
factor upf. 

These states are exact eigenstates of SU(4) invari- 
ant Hamiltonians at the integer partial filling factors 
I'N — 0,1,2,3, but are they the ground states? Quan- 
tum Hall ferromagnetism theories are based on the no- 
tion that they normally are. Whether or not this idea 
is correct likely depends to some degree on the charac- 
ter of the Hamiltonian, with Hamiltonians that are more 
strongly repulsive at short distances most likely to sup- 
port these ground states. In the case of SU(2) (spin flavor 
only) , it seems clear from finite-size exact diagonalization 
calculation's!! that they are ground states for small n 
Landau levels. Exact diagonalization calculations for the 
graphene SU(4) case also support the quantum Hall fer- 
romagnetism ansatz for small N in graphene. Although 
it is not possible to prove rigorously that the integer i^n 
ground states have this form, even when Landau level 
mixing is neglected, the case for quantum Hall ferromag- 
netism at small |7V| appears to be quite strong. Its ne- 
glect of Landau-level mixing effects also appears to be 
soundly based although this approximation does not be- 
come better, as it does in ordinary 2DEGs, as the field 
strength is increased. 

When the quantum Hall ferromagnetism ansatz is cor- 
rect, an electron that is added to an integer partial fill- 
ing factor state will have to occupy a flavor Landau level 
that is empty, losing exchange energy in the process. The 
charge gap is due to this loss of exchange energy. This in- 
compressibility at a fixed filling factor leads to an integer 
quantized Hall conductivity and to vanishing longitudi- 
nal conductivity, in accord with well-established quan- 
tum Hall phenomenology. 

Since SU(4) is a continuous symmetry, when it is bro- 
ken spontaneously the systems must support gapless col- 
lective mode s. Th ese neutral, spin-wave like modes have 
been studied,'^212S] and their full spectra have been deter- 
mined exactly.'^ Because the symmetry group is larger 
than in the SU(2) case,'^ the number of spin-wave-like 
modes is larger than one. The number of modes is given 
by J^Ar(4 — vn) — four modes for = 2 and three 
modes vp^ — 1, S.'^Sl Perhaps more interesting is the 
appearance of topological solitons in the order parame- 
ter called Skyrmions as low-energy charged excitations 
of the systems.'^^'^ Again due to the larger symmetry, 
Skyrmions come in more species in the SU(4) quantum 
Hall ferromagnetj2l! case compared t o the SU(2) quan- 
tum Hall ferromagnets studied earlier .I^IEIl 

There is an important quantitative difference com- 
pared to the ordinary 2DEG case in the spatial struc- 
ture of electron-electron interactions that is due to the 



Dirac nature of graphene's charge carriers. In a mag- 
netic field, the two sublattice components of graphene's 
Dirac-model eigenspinors have different orbital Landau 
level wave function (gee also Section IV A). As a 
result, the form factor which determines the effective 
electron-electron interaction in a given Landau level is 
quite different for Dirac fermions.'^SESl One of the con- 
sequences of this is that Skyrmions are the lowest en- 
ergy charged excitations not only in the lowest Landau 
level (which is the only stable Skyrmion case for non- 
relativistic el ectron jSH) ^ but also in some of the higher 
Landau leveld^IlM j^ore specifically they are lowest en- 
ergy charge excitations for |A^| < 3.^— These high Landau 
level skyrmions have been observed in recent numerical 
work.l^ 

Because states in different valleys are certain to be on 
different sublattices in the iV = case, valley and sublat- 
tice (layer in the bilayer case) labels are equivalent when 
Landau level mixing is neglected. This situation changes 
when mixing becomes important. The role of Landau 
level mixing can be increased in the single-layer case by 
increasing the strength of carrier-carrier interactions, i.e 
by decreasing dielectric screening. (This contrasts with 
the bilayer case in which Landau level mixing can be en- 
hanced simply by decreasing the magnetic field strength.) 
When mixing is important there is a qualitative difference 
between sublattice and valley labels because sublattice is 
not a good quantum number. When the valence Lan- 
dau level breaks symmetry by establishing a condensate 
which favors sublattice polarization, the wavefunctions of 
Landau levels far from the Fermi surface are altered in a 
way which reinforces the condensate. In the language of 
the 2D Dirac equation, a sublattice polarization term is a 
mass term, as discussed further below. Depending on the 
relative signs of the masses for different spins and valleys, 
sublattice symmetry breaking leads to various kindsPof 
density wave states and to gaps at a variety of integer 
filling factors that persist to weak fields. Quantum Hall 
ferromagnet states in which the broken symmetry con- 
densate produces a mean field term which favors occu- 
pation of one sublattice over another and plays the role 
of this mass in the 2D Dirac equation. For total Landau 
level filling factor v = this subset of quantum Hall fer- 
romagnet broken symmetry states is related to the idea 
of magnetic catalysis. Indeed there is a subset of quan- 
tum Hall ferromagnet states which, when stable, tend to 
remain robust at weaker magnetic fields. There is some 
evidence that the v = Q state in both single-layer and 
bilayer graphene may have the type of order anticipated 
by the magnetic catalysis scenario. 



2. Explicit SU(4) Breaking Interactions 

A number of small perturbations that explicitly break 
SU(4) symmetry are present in graphene. They can fix 
the orientation of the SU(4) order parameter, i.e. select 
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specific preferred flavor orbitals to occupy at a particular 
partial filling factor, and will also open small gap(s) in 
some of the the collective mode dispersions. In the fol- 
lowing we discuss some of these perturbations. 

Zeeman Splitting: — This is probably the most obvious 
explicit symmetry breaking interaction. It lifts the spin 
degeneracy without affecting the valley degeneracy; for- 
mally it reduces the SU(4) symmetry to a (valley) SU(2) 
(8) U(l) symmetry. Its influence is particularly important 
when i/N = 2, in which case it uniquely selects a single 
ground state which is a valley singlet but has spins in 
the valence Landau level fully polarized along the field 
direction. In the absence of interactions the charge gap 
associated with this state is precisely the Zeeman split- 
ting, and in some theories (to be discussed later) this is 
thought to be the sole source of gap ai v = ±4, when 
the = 1 Landau level is half-filled. Quantitatively, the 
Zeeman splitting is of order 10 K at 10 Tesla, compared 
to the electron-electron interaction scale e^/e^ ~ 500if 
{I — \fhcfeB is the magnetic length) ; it thus qualifies as 
a small symmetry breaking perturbation. 

Lattice Effects: — SU(4) symmetry is exact in the con- 
tinuum limit a ^ £, where a is the lattice spacing. 
Lattice e ffects break SU(2) symmetry within the valley 
subspace.'SSEZlThe strength of this symmetry breaking is 
proportional to a/£. For the lowest Landau level (TV — 0), 
lattice effects introduce an Ising anisotropy in the val- 
ley SU(2) spacCjl^^ while for \ N\ > an easy-plane (or 
XY) anisotropy is introduced.'^SEH Lattice Hartree-Fock 
calculation*^ point to one important role of the lattice. 
In a continuum model electrons do not change their valley 
index when they scatter off another electron. As a conse- 
quence exchange interactions occur only within valleys. 
Although the neglected inter-layer exchange interactions 
are weak, they do favor sublattice polarization states in 
which the sense of layer polarization is the same for both 
valleys. This inter-valley exchange effect can be impor- 
tant in determining the character oiv = Q quantum Hall 
states as we discuss further below. 

Disorder: — Although disorder potentials that are 
smooth on atomic length scales do not break any symme- 
try, they do lift Landau level degeneracy. Disorder favors 
SU(4) singlet ground states because more electrons can 
be placed where the disorder potential is lowest when 
the four valley degrees of freedom are exploited symmet- 
rically. In Ref. [25l in which the notion of QHF was first 
introduced in the graphene context, a Stoner-like crite- 
rion was derived for the competition between interactions 
and disorder in quantum Hal ferromagnetism. The situ- 
ation was also analyzed in Ref. [27l in which the authors 
came to the conclusion that the samples available at that 
time might not be clean enough for QHF, and the systems 
under study at that time, and perhaps now, might actu- 
ally be SU(4) paramagnets. They then argued that the 
combination of Zeeman and symmetry breaking interac- 
tions can give rise to new QH states in the paramagnetic 
regime. 

On the other hand Abanin et al.^^ pointed out an inter- 



esting symmetry breaking effect of the disorder potential, 
which may favor ordering: a random potential produced 
by strain in the graphene lattice locally breaks the val- 
ley degeneracy in a random manner; it thus acts like a 
random field along the z direction in the valley subspace. 
This could force the valley SU(2) order parameter into 
the XY plane, and allows for algebraic long-range order 
at finite temperature and a Kosterlitz-Thouless transi- 
tion. 



3. Magnetic Catalysis 

One set of theorie^^HllIllEl] related to interaction 
induced gaps in graphene's quantum Hall regime fo- 
cuses on interaction induced mass terms in the two- 
dimensional Dirac equation that describes graphene elec- 
trons. Interaction-induced masses in the absence of 
a magnetic field are partly analogous to interaction- 
induced masses in elementary particle physics. Mag- 
netic catalysis refers to the reduction in the interac- 
tion strength necessary for spontaneous mass that occurs 
when a magnetic field is present. 

Since a Dirac equation mass corresponds to a CTz term 
in the sublattice-pseudospin mean-field Hamiltonian, it 
is always generated in states that have A-B sublattice 
polarization within each valley. The broken symmetry 
states anticipated in magnetic catalysis theories are ones 
in which the same sense of sublattice polarization oc- 
curs in each valley. From a physical point of view these 
are therefore density-wave states in which the honeycomb 
lattice symmetry is reduced. Spontaneous mass genera- 
tion in monolayer graphene was in fact discussed as a 
theoretical possibility even before the experimental real- 
ization of the material.!^ The excitonic gap imagined in 
this early work was attributed to an instability towards 
electron-hole pairing in the presence of a constant mag- 
netic field in (2 -I- 1)D-QED. Magnetic catalysis theories 
are distinguished phenomenologically from QHF theories 
by their prediction that spontaneously generated gaps 
should appear in the N = Landau level only. Induced 
mass type order opens gaps between the conduction and 
valence bands, instead of opening them within all Lan- 
dau levels. When Zeeman splitting is included, these 
theories predict plateaus in the Hall conductivity axy at 
1^ = 0, ±1, and I' = ±2k for fc = 1, 2, 3, • • • . In Ref. El 
it was argued that the ordering that generates the mass 
could be either a charge density wave (CDW) or a spin 
density wave (SDW), depending on the relative strength 
of on-site and nearest neighbor electron-electron repul- 
sion. SDW states have the opposite sense of sublattice 
polarization for opposite spins. 

Earlier experimental resultJ^^ appeared to be consis- 
tent with the magnetic catalysis scenario since the first 
odd integer QH states to be observed appeared at = ±1 
in the = Landau level. (Assuming that spins are 
maximally polarized because of the Zeeman field, gaps 
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at v = ±1 require that valley or sublattice symmetry be 
broken, consistent with the magnetic catalysis scenario.) 
However the recent observation of integer QH states at 
I' = 3 in suspended graphene samples and at all integer 
filling factors for graphene on h-BN substrate seem to 
indicate the QHF scenario better describes the physics 
which takes place in graphene sheets. In fact magnetic 
catalysis should be viewed as a special case of quantum 
Hall ferromagnetism in which the broken symmetry state 
is a charge-density-wave or spin-density-wave state. 

CDW and SDW states are likely strong-interaction- 
limit broken symmetry states for honeycomb lattice mod- 
els in the absence of a magnetic field. If these broken 
symmetries were incipient at B = 0, a weak magnetic 
field would in all likelihood push them over the edge — 
hence magnetic catalysis. If on the other hand B = 
graphene systems are not close to instabilities, there is no 
reason to expect that CDW and SDW instabilities would 
dominate in the presence of a magnetic field. The exper- 
imental observation of QHE's at nearly all integer filling 
factors suggests that it is Landau-level formation and 
interactions within groups of degenerate Landau levels, 
the quantum Hall ferromagnetism scenario, that plays 
the essential role in establishing these interaction-induced 
gaps. If so, this piece of graphene phenomenology sug- 
gests that the B — graphene system is not especially 
close to an interaction-induced instability. However in bi- 
layer graphene it appears that the story is different. Some 
quantum Hall effects appear to persist to very weak mag- 
netic fields^^^^° and sometimes to B = 0, with interac- 
tions playing an essential role. It appears likely that neu- 
tral bilayer graphene band structure is unstable at B = 
against weak interactions, as predicted theoreticallji^^^'^ 
(see section V. B). 



4- = Quantum Hall State Transport Properties 

The v — QH state in graphene has garnered par- 
ticular interest since it does not have an analogue in 
semiconducting 2DEGs with parabolic dispersion. The 
new physics is due to the relativistic Dirac dispersion for 
quasiparticles in graphene which lead to conduction and 
valence bands (see Section HI A) that are not separated 
by a gap. Initial experimental studies of this state in- 
dicated that the transport phenomenology on the v — 
QH state is different from that of the other integers. This 
unusual behavior is evident both in p^x and Pxy, the diag- 
onal and off-diagonal resistivity respectively. The exper- 
imentally measured Hall resistivity does not display ei- 
ther a clear quantized plateau or a zero resistance statef^ 
however a step- like feature was evident in the off-diagonal 
conductivity at a^y — 0. The longitudinal resistivity did 
not exhibit a deep minimum normally characteristic of 
the conventional QHE. Later experiments on the ly — 
state reported a steeply divergent Dirac point resistance 
at low temperatures.^^ The authors argued that this in- 
sulating behavior followed Kosterlitz-Thouless (KT) type 
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FIG. 6. Schematic representation of the energy dispersions 
of edge states of (a) valley-split {Av > As) vs (b) spin-split 
(As > Ay) V = QH state. In the valley-split scenario 
(a) there are no gapless edge states and the system resem- 
bles an insulator. Whereas the spin-split scenario supports 
counterpropagating spin-polarized edge states with the sys- 
tem identical to the quantum spin Hall state. 



scaling with respect to the magnetic field and, was inter- 
preted as a quantum-KT phase transition by assuming 
a valley split XY state. This field induced insulating 
behavior has also been observe d in s uspended graphene 
samples in a two terminal setup .^^SEl] 

The unconventional transport phenomenology of the 
v — QH state has drawn interest from a number of 
theorists. Most work can be divided into two camps, 
one of which concentrates on the edge state structure 
whereas the other focuses on the bulk properties of the 
V = state. Due to the inherent particle-hole symmetry 
of the Dirac equation, the edge states in graphene®^ are 
remarkably different from those in conventional 2DGEs. 
Near the physical edge of the sample, the particle-like 
branches disperse upwards in energy whereas the hole- 
like branches disperse downwards as indicated in Fig. [5] 
For the N = Landau level particle-hole symmetry dic- 
tates that each spin has one electron-like and one hole- 
like edge state branch. In the presence of interaction 
induced gaps the edge state structure would be different 
for a valley- or spin-split ly — QH state. The valley-split 
state can be truly insulating when the Fermi energy lies 
in the gap since there would be no gapless edge states as 
shown in Fig. [6^. On the other hand a spin-split state 
would sup port c ounter-propagating edge states indicated 
in Fig. [6|d,IS^^ a scenario identical to that of quantum 
spin Hall systems.^^ The unusual transport phenomenol- 
ogy of the v = QH state was initially explained as- 
suming an interaction enhanced spin-split gap.^^ How- 
ever the observation of the field induced insulating fea- 
ture cast some doubts on the counter-propagating edge 
state picture. Shimshoni et. al^ proposed that scatter- 
ing between the counterpropagating edge modes due to 
the presence of localized moments can lead to a metal- 
insulator quantum phase transition, saving the spin- 
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polarized scenario. In this picture the crossover from 
metaUic to insulating behavior follows a quantum KT 
type scaling in the magnetic field. 

The observation of the a divergent resistance in Ref 
[76l along with the KT scaling raised questions about the 
completeness of the edge state picture. Based on the 
experimental results it was argued the the v — QH state 
must exhibit XY valley-pseudospin order, where phase 
coherence is established between valley K and K'.l^ It 
is well known that this state supports vortex and anti- 
vortex topological excitations. It was argued that in this 
case unbinding of the vortex-anitvortex pairs triggers a 
bulk KT transition from the pseudospin XY quasi-long 
range ordered phase to the disordered phase at a critical 
magnetic field Be- 
Most theoretical studies of the v — QH state have 
been performed based on the Dirac-equation continuum 
models. These recent experiments raise the suspicion 
that the ground state for the ly — QH state could be in- 
fluenced by physics that is outside of the Dirac-equation 
continuum model. Calculations based on 7r-band-tight- 
binding model including electron-electron interactions 
performed by one of us indicate that the ground state 
for graphene sheets on substrates is most likely a spin- 
density wave, and that a charge-density-wave is possi- 
ble for suspended samples.'^ These states are favored by 
inter-valley exchange interactions. 

It was pointed out in Ref. HO! that one way to distin- 
guish between edge and bulk mechanisms of this appar- 
ent metal insulator transition is to study samples with 
Corbino geometry, in which edge states play no role. 
Thus the same transition should be observed in Corbino 
samples as magnetic field increases, if the transition oc- 
curs in the bulk. On the other hand if the transition is 
only an edge phenomenon, then Corbino samples should 
only exhibit insulating behavior. Ref. 80 further dis- 
cusses the possible role that disorder plays in a bulk 
transition: disorder induces inhomogeneous charge dis- 
tribution, which results in puddles with different densi- 
ties and thus different filling factors in the presence of 
magnetic field; edge states separating different puddles 
may play an important role in bulk transport. Increas- 
ing magnetic field suppresses such puddle by increasing 
spacing between different Landau levels; this provides a 
possible bulk mechanism of the transition. This scenario 
appears consistent with experimentally observed charge- 
inhomogeneity in the presence of magnetic field.'^ 



B. Interaction induced integer QH states in bilayer 
graphene 

The concept of QHF is quite general and applies 
whenever Landau levels with some internal flavor indices 
are degenerate. The N ^ Landau levels of bilayer 
graphene have the same spin/vallcy four-fold degeneracy 
as single layer graphene and can be described by an 
approximate SU(4) invariant Hamiltonian. The ground 



states corresponding to the interaction induced plateaus 
are then captured by the variational wavefunction (21 1. 



Generally the physics associated with these states should 
be similar to that of the higher Landau levels in single 
layer graphene, with the same number of collective 
modes and the presence of Skyrmions. However there 
will be some important quantitative differences related 
to bilayer graphene's orbital structure (see section IV 
A. Eq. 18 1. This is because the spatial structure of 
the electron-electron interactions is modified due to 
the J = 2 chiral nature of bilayer graphene's charge 
carriers.l^SI In particular Skyrmions are the lowest energy 
charge excitation for \N\ < 2 in bilayer graphene. 
One important difference between single and bilayer 
graphene is that in the bilayer case it is easy to generate 
a mass terms (an energy difference between sublattices 
Ai and A2) externally, simply by applying an electric 
field between the layers.- In the N ^ Landau levels 
this should allow the study of transitions between the 
different broken symmetry states function of an 

external control parameter. 

The most interesting aspects of bilayer graphene 
are associated with the interaction driven states that 
form between v = —4 and v = A QH plateaus ,1^2121] 
many of which do not have a counterpart in either 
graphene or other semiconducting QH systems. The 
iV = Landau level manifold in bilayer graphene 
has an extra pseudospin corresponding to the n = 
and n = 1 orbital degree of freedom.!^ The ordering 
and fluctuations of this pseudospin degree of freedom 
is unique to bilayer graphene and leads to unusual 
collective modes,^ anomalous exciton condensation,!^ 
and spiral ordering.!^ It is this new aspect associated 
with the physics of bilayer graphene that will be the 
focus of the rest of this section. 

The iV = Landau level of neutral bilayer graphene is 
the direct product of three S* = 1/2 doublets: real spin, 
which layer pseudospin as in semiconducting double 
well quantum Hall systems (DWQHS) and the Landau 
level pseudospin n = and n — I which is responsible 
for the new physics. When the magnetic field is strong 
enough or disorder is weak enough, interactions lead to 
quantum Hall plateaus at the octet's intermediate inte- 
ger filling fact ors b oth for suspended bilayer^^ and on 
SiO substratesl ^^l^^l For neutral bilayers the octet filling, 
proceeding in integer increments starting from the filling 

— —4, can be characterized by Hund's rule behavior 
of the pseudospin: first maximize spin-polarization, 
then maximize layer-pseudospin to the greatest extent 
possible and finally maximize the Landau level pseu- 
dospin polarization to the extent allowed by the first 
two rules.!^ For balanced bilayers interactions favor 
states with spontaneous coherence between the valleys, 
and hence layers, since these broken symmetry states 
do not require charge transfer between the layers.!^ 
The Hund's rule behavior indicates that all odd integer 
filling factors between v = —4 and v — 4 have Landau 
level pseudospin polarized. This preference to polarize 
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the n = Landau level pseudospin is because of a larger 
gain in exchange energy over the n = 1 Landau level 
pseudospin and is related to the wavefunction profile of 
the n = and n — 1 Landau levelP^l 

The octet state in bilayer graphene supports 
i'Ar(8 — vn) collective modes with i^n defined simi- 
lar to that in section V A. The orbital fluctuations 
associated with the Landau level-pseudospin are dom- 
inant at odd fillings of the octet. Exchange energy 
consideration at even integer fillings require that the 
Landau level-pseudospin degree of freedom gets frozen 
out, hence the interaction driven states are more similar 
to the semiconducting DWQHS. One particularly inter- 
esting and unusual collective mode associated with the 
balanced bilayer at the odd fillings is the intra-Landau 
level mode.**^ This collective mode exhibits a gapless 
'^q ^ Q^^^ dependence in the long wavelength limit, a 
behavior that is due to the presence of long range dipolar 
interactions associated with the Landau level pseudspin 
fluctuations. Fluctuating Landau level pseudospinors 
are a linear combination of n = Landau level orbitals 
(even with respect to their cyclotron orbit center) and 
n = 1 orbitals (odd with respect the orbit center), and 
therefore carry an electric dipole proportional to the 
in-plane component of their pseudospin. Since the dipo- 
lar interactions are long ranged they play a dominant 
role at long wavelengths leading to the non-analytic 
long wavelength behavior.l^ The application of Kohn's 
theorem within the projected lowest Landau level leads 
to a gapless mode due to a vanishing contribution from 
electron-electron interactions, hence this intra-Landau 
level mode would acquire a gap in the presence of a 
single particle induced Landau level-splitting. It was 
proposed in Ref. that this mode could be observed 
in cyclotron resonance signal by applying an external 
potential difference between the layers. The external 
potential leads to Landau level-splitting which is positive 
(n = 1 above n = orbital) in one valley and negative 
(n = 1 below n = orbital) in the other.ES 

The presence of the Landau level-pseudospin also 
has important consequences on the properties of the 
exciton condensates at odd fillings v = —3, 1. The 
superfluid density vanishes at these filling factors 
leading to a finite-temperature fluctuation-induced first 
order isotropic-smectic phase transition when the layer 
densities are not balanced. This vanishing superfluid 
stiffness leads to a quadratic rather than expected 
linear'^ phonon collective mode. The transition to the 
smectic state is a consequence of the negative Landau 
level-pseudospin gap in one of the valleys along with 
the zero superfluid stiffness; combined this leads to a 
long-wavelength instability of the Brazovskii-type.'^^'^ 
At other odd filling factor u = —1,3, due to a neg- 
ative Landau level-gajPSl^ the single particle Landau 
level splitting competes with the interaction induced 
Landau level exchange splitting which prefers n = 
Landau level polarization. This leads to a sequence of 
transitions from i) an inter-layer coherent state, ii) a 



mixed state with inter-orbital and inter-layer coherence 
to iii) an inter-orbital coherence state spiral Landau 
level-pseudospin ordering.!^ The spiral state is due to 
the presence of Doyzaloshinskii-Moriya terms which is 
related to the inversion symmetry breaking of the QH 
ferromagnet. The presence of orbital coherence leads to 
a finite density of in-plane electric dipoles^ which can 
be manipulated by in-plane electric fields. 

Another interesting aspect associated with the broken 
symmetry = octet states in bilayer graphene is the 
possibility of new and interesting topological excitations 
not expected in graphene and semiconducting 2DEGs. 
It has been proposed that charge-2e Skyrmions are the 
low energy charged excitation near the filling factors 
I' = —2, 2.*^ As expected these charge-2e Skyrmions 
have textures in both n = and n = 1 Landau level 
orbitals. Interlayer charge-2e meron and Skyrmions 
crystals with checkerboard patterns for balanced and 
unbalanced layers have subsequently been identified in 
numerical calculations.!^ These charge-2e objects also 
condense in a rich variety of crystal states at odd filling 
factors where these pseudoskyrmions are associated with 
the Landau level-orbital degree of freedom. The orbital 
Skyrmions exhibit an unusual vorticity and charge rela- 
tionship exhibiting textures of in-plane electric dipoles 
and can be seen as an analogue of spin Skyrmions which 
carry a magnetic textures. These texture would couple 
to an in-plane electric field and the modulation of the 
electronic density in the crystalline phases should be 
experimentally accessible through a scanning tunneling 
microscope measurement of their local density of states. 

Just like monolayer graphene v = state in bilayer 
graphene has drawn special attention. In the quantum 
Hall regime the ordering see ms to indicate either spin 
or valley polarization, or both.l^Slill xhe interest in the 
1^ = QH state of bilayer graphene is different from that 
in graphene though, as the parabolic dispersion of bilayer 
graphene makes it susceptible to interaction-induced 
symmetry breaking even at zero magnetic field.'^^lt^ 
Incompressibilities for suspended bilayer graphene were 
seen to survive even in the B — limit, where it 
was concluded that the incompressibility at the charge 
neutrality was consistent with either spontaneously 
broken time-reversal symmetry or spontaneously bro- 
ken rotational symmetry. The broken time reversal 
symmetry leads to an anomalous quantum Hall state 
where electrons in one layer occupy valley K whereas 
electrons in the other layer occupy valley K'.'^ This 
state supports topologically protected edge modes and 
would exhibit a finite conductance. The time-reversal 
symmetric state is predicted to exhibit nematic order 
where the parabolic dispersion is expected to split into 
two Dirac cones. This would result in a lowering of the 
density of state s at th e charge neutrality point 

Some author^^^Ml have concluded that the B = 
broken symmetry states have spin and valley dependent 
sublattice polarization. If this is true the states are 
associated with interesting momentum space Berry 
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curvaturej^ anomalous Hall, and orbital magnetism 
effects. Among the states in this class with no overall 
layer polarization, ar^^l SDW states, usually called LAF 
states (layer antiferromagnet) in the bilayer case, QSH 
(quantum spin Hall) and QAH^'' (quantum anomalous 
Hall states). These three states have quantized Hall 
conductances with the values 0, 0, and ±4e^//i for 
the LAF, QSH, and QAH states respectively. In each 
case the gaps associated with order evolve smoothly 
into quantum Hall gaps with filling factors v = 0, 0, 4 
respectively. In a magnetic field the QAH state evolves 
into a bilayer state that is similar to the non-interacting 
I' = A state. The QSH state evolves in field into a 
spin-polarized state with total filling factor = but 
a filling factor difference of 4 between majority and 
minority spins. The LAF state evolves into a state 
which appears at v = and has no spin-polarization. 
Given these properties, studies of the quantum Hall 
effect at weak magnetic fields can shed light on the 
character of the B = broken symmetry states. For 
example if the QAH state is the B — ground statd^ 
the V = A quantum Hall effect should persist to S = 0, 
whereas if the LAF state where the ground state it is the 
1/ = quantum Hall state that would persist. Current 
experiments do not yet provide consistent results on 
these points, presumably because of small differences in 
disorder potential. The possibility of observing quite 
different transport properties on quite similar samples 
suggest that the various potential ordered states of 
bilayers compete closely. 



VI. FRACTIONAL HALL EFFECT IN 
GRAPHENE 

As explained earlier electron-electron interactions in 
semiconducting 2DEGs can give rise to incompressibili- 
ties at fractional values of the Landau level filling factor 
V. Interactions within a partially filled Landau level lead 
to the formation of strongly correlated states that exhibit 
fractionally charged excitations and fractional statistics"* 
which in some cases are proposed to be non-Abelian.l^ 
Graphene is very similar to its semiconducting 2DEG 
cousin; in fact there are more similarities than differences 
as far as FQHE are concerned. The two dimensionality. 
Landau level structure and the appearance of integer QH 
effects due to both Landau quantization and interaction 
induced gaps are good benchmarks indicating that in- 
compressibilities at fractional fillings are likely to exist. 
While similarities do exist there are also differences, due 
to the added valley degeneracy and the possibility of tun- 
ing Coulomb interactions strength by substrate manipu- 
lation. As a result of these features FQHE in graphene 
is expected to be richer. Specifically the four-fold de- 
generacy in graphene can lead to spontaneously broken 
SU{A) symmetric FQH states with an ordering of inter- 
nal indices that lead to a multitude of FQH states not 



observed or anticipated in semiconducting 2DEGs.l^ 

FQH physics is expected to dominate when electron- 
electron interactions are stronger than disorder poten- 
tials. At the time of this writing it still appears to be 
true that this ratio has a larger value in some semicon- 
ducting quantum well systems than in the cleanest avail- 
able graphene samples. Current theories of FQHE in 
semiconducting systems rely heavily on the notion of the 
projection of interactions onto a single Landau level (i.e 
in general cases Landau level mixing is neglected). The 
appropriateness of the Landau level projection is not triv- 
ially obvious for the case of graphene. In semiconducting 
2DEGs this is formally achieved in the limit of a large 
magnetic field B oo. This is justified because inter- 
action gaps T4„e ^ /els scale as \/B while the single 
particle Landau level gap hiOc scales as -B. So it can 
be argued there that for large magnetic fields interac- 
tions between electrons in the lowest partially-filled Lan- 
dau level cannot induce transitions to the higher Lan- 
dau levels. The projection of the interactions onto the 
lowest partially-filled Landau level is not as clearly justi- 
fied in graphene because the single particle band gaps in 
graphene also scale as ~ \/B, which is the same as the 
interaction strength. However, the validity of Landau 
level projection in theories of the graphene FQHE can 
be justified by considering the large dielectric constant 
limit in the surrounding material. Most FQHE theory 
for graphene has nevertheless assumed that Landau level 
mixing is an inessential complication that can be ignored. 

On the experimental side, the focus in graphene has 
shifted from IQHE to FQHE since higher mobilities have 
been achieved both by suspending samples and using h- 
BN as a substrate. The FQHE was initially observed 
in suspended samples via two terminal transport mea- 
surements which revealed clear plateaus corresponding 
to the V = 2 — 5/3 = 1/3 FQH state developing as a 
function of the magnetic field.'^^*^ More recently a se- 
quence of FQH states have been observed in high mobil- 
ity multi-terminal graphene devices fabricated on a h-BN 
substrate.'^ The observed fractions are all thirds (i.e. of 
the form v = fc/3) and belong to the = and A^ = 1 
graphene Landau levels. In the A^ = Landau level, the 
observed fractions indicate multicomponent FQH states 
with spontaneously broken SU (4) symmetry, and possi- 
bly full symmetry broken states in A'^ = 1 Landau level. 
It is interesting to note that none of the fractions cur- 
rently observed in the A^ = Landau level of graphene 
have an analogue in the semiconducting GaAs; for ex- 
ample none of the observed fractions correspond to the 
celebrated Laughlin state at = 1/3. The observed 
fractions are most likely multicomponent generalization 
of the Laughlin state. Variational wavefunctions with an 
internal SU{2) symmetry were introduced by Halperin 
by generalization of the Laughlin state to spin-1/2 par- 
ticles. In the following we review the essential properties 
of the Laughlin state and provide the SU(4:) generaliza- 
tion of these states which are most likely related to the 
FQH states observed recently. 
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The FQHE is characterized by the condensation of 
electrons into special highly correlated states that min- 
imize the Coulomb energy. Interactions play a crucial 
role in determining the nature of these correlated states. 
A fractionally filled Landau level of noninteracting elec- 
trons does not define a unique many-body state since 
electrons can be placed into any of the degenerate or- 
bitals in an exponentially large number of ways subject 
only to the Pauli principle. Interactions lift this degen- 
eracy by selecting a unique ground state from the vast 
space of degenerate states, in which electrons avoid each 
other as much as possible. At some specific filling factor 
i^N, the energy change upon adding one electron differs 
from the energy change for removing one electron by a 
finite charge gap. The resulting fiuid is incompressible 
and often characterized by unusual topological order.^^l 
For fractional fillings of the form i^n — l/m the essence 
of this phenomenon is captured by Laughlin's wavefunc- 
tion, 



.]=n(^.-^.rexp(-i^|z,p), 



(22) 



where Zi — Xi + iyi denotes the two-dimensional positions 
of the electron in terms of a complex number, and we have 
set the magnetic length = 1. Since this wavefunction 
lives in the lowest Landau level (for = 0, correspond- 
ing wave functions exist for any N) the parameter m 
must be an integer to satisfy wavefunction analyticity.l^ 
Furthermore to satisfy the wavefunction antisymmetry 
requirement for fermions, m must be odd. The Laugh- 
lin wavefunction has good correlations built in it since 
each electrons sees m-zeros at the positions of all the 
other electrons. This helps minimize the Coulomb in- 
teractions uniformly among all the electrons. The prob- 
ability of two electrons approaching each other, which 
is given by the amplitude of the wavefunction, vanishes 



rapidly: p{r,-rj) ^ \ri-r. 



1 2m 



Note that no pair of par- 



ticles ever has relative angular momentum less than to, 
and therefore the Laughlin state is a zero-energy exact 
eigenstates of model Hamiltonians initially introduced by 
Haldane in which electron only interact if they have small 
relative angular momentum. In the following we briefiy 
review this point of view on FQH physics for semicon- 
ducting 2DGEGs and apply these ideas to describe wave- 
functions for SU (4) FQHE in graphene. 

In semiconducting 2DEGs the primary features of FQH 

(n) 

states are captured by Haldane pseudopotentials Vm 
which represent the interaction energy of two electrons 
with relative angular momentum to. Haldane's pseu- 
dopotentials are essential for exact diagonalization cal- 
culation for small sized systems and can give insight into 
possible FQH states by looking at overlaps of trial wave- 
functions with exact numerical solutions. Furthermore 
this procedure introduces model Hamiltonians for which 
exact eigenstates which exhibit incompressibility at frac- 
tional filling are available. An example of this is the hard 
core model for which the zero energy eigenstate is the cel- 
ebrated Laughlin state. Any given interaction y(|ri — rj|) 



can be projected onto relative angular momentum basis 
|m') where to', 



i<j ra' 



where 



n'\V\m') 
(to'|to') 



(23) 



(24) 



are known as the Haldane pseudopotentials and P™ is 
the projection operator which selects particles with rel- 
ative angular momentum m! . The Laughlin state is the 
exact zero energy eigenstate of the "hard-core" potential 
with Vm' = for to' > to. This simply follows from the 
fact that \il'm) = for any m' < m since every pair 
has relative angular momentum of at least to. These 
models allow for a simple generalization to the case of 
multiconiponcnt wavefunctions. 

By mapping Hamiltoninans projected onto different 
isolated Landau levels, Haldane pseudpotenetials can be 
calculated for any value of N. In the A^ = Landau 
level in graphene the wavefunctions are localized on sub- 
lattice A{B) for valley K(K') with orbital index n — 0, 
whereas for the N graphene Landau level the wave- 
function is a coherent superposition of n and n — 1 or- 
bital wavefunction of semiconducting 2DEGs localized of 
different sublattices. It follows tha ti^^'^^' the A^ = Lan- 
dau level graphene pseudopotentials are identical to the 
lowest Landau level pseudopotentials in semiconducting 
2DEGs, however the A^ ^ pseudopotentials are differ- 
ent and given by. 



(27r) 



[L|jv|(gV2)+i|jv|-i(gV2)] 



(25) 

where Ln{x) is the n*'* order Laguerre polynomial and 
Vq is the Fourier transform of the interaction potential 
with Vq = 27re^/(eg) for Coulomb interaction. The rel- 
ative strength of the Haldane pseudopotential /V^ 
normally indicates the possibility of finding strongly cor- 
related FQH states. Numerical investigations have shown 
that graphene should support FQH-like correlated states 
in the lowest Landau level.'^S 

Motivated by Haldane pseudoptential considerations 
one can write a class of model Hamiltonians with exact 
zero-energy eigenstates. Relevant to graphene are the 
SU{4) generalization of the Halperin-Laughlin wavefunc- 
tionpl 



= n n 

a — 1 ia <ja 



n n (4 

a<P lad 13 



P \nal3 



(26) 

where we have omitted the ubiquitous Gaussian factor. 
In the above equation — x"^ + iyj^ denotes the po- 
sition of the i*'^ particle with a flavor index a. It is easy 
to write the relevant model Hamiltonians for which these 
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wavefunctions are the zero-energy eigenstates, we leave 
this as an exercise for the reader. Apart from the omit- 
ted exponential factor, the wavefunction must be a ho- 
mogeneous polynomial of degree in z"^ where 
is the degeneracy of of the a Landau level. Further- 
more the polynomial must be antisymmetric under ex- 
change of and Zj^ for the same index a, which re- 
strict rUaS to take odd integer values. Using the fact 
that every extended single particle orbital in the lowest 
Landau level involves a polynomial with a fixed density 
of zeros given by B/^q (where 3>o = hc/e is the flux 
quantum) one can express the density pa of each a com- 
ponent i?/$o — PaiTia + ^a^ii Panafj. The individual 
filling factors can be calculated using 



(27) 



by defining a symmetric K matrix with the diagonal 
elements kaa = "m-a and the off-diagonal components 
kap = Tiap- The total filling factor corresponding to K 
is given by vn = Vi. Eq(27l only makes sense if the 
matrix K is invertible. One can classify the spin and 
valley polarization of these states based on the rank of 
the corresponding K matrix and the filling factors. For 
example since a rank 4 matrix is invertible the individ- 
ual spin and valley flavor filling factors Vi are uniquely 
determined, therefore the spin polarization in the z- 
direction is just given as a 5^ = X^o-l^o-.T ~ ^a-i) 
with a similar expression for the valley polarization.^^ 
The single component v = \ jm Laughlin like states cor- 
respond to K matrices of rank 1 with all the exponents 
ma = JT-ce,^ = JTi, such 'A State would be a fully polarized 
SU(4)-ferromagnetic state. The consideration of spin and 
valley ordering is important in determining the ground 
state for a filling factor as quantum Hall ferromagnetic 
physics prefers states that exhibit ferromagnetic order- 
ing. The if-matrix can be viewed as a SU(4) generaliza- 
tion of Wen's K-matrix formulation of the FQH effect. 
This ii'-matrix also encodes the quasiparticle statistics 
and charge for a given state. 

From the experimental results on h-BN substrate it 
seems that these SU(4) states describe the observed 
fractional plateaus within the lowest Landau level in 
graphene. It seems that in the lowest Landau level the ex- 
plicit symmetry breaking terms including valley and Zee- 
man splitting are sufficiently small that the Coulomb in- 
teractions intermix the four branches, allowing for mixed 
valley-spin FQH states.l^ One can argue this based on 
the pattern of the particle- hole symmetry of the observed 
FQH states; for example particle-hole conjugate FQH 
state must share similar gaps. The observation of the 
fractional states v = 1/3, 2/3 and 4/3 and their conjugate 
fractions seem to indicate that the external symmetry 
breaking terms are sufficiently small and an approximate 
SU(4) symmetry is preserved. The FQH states observed 
in the = 1 Landau level do not follow this patterns 
indicating that all degeneracies are lifted either due to 
explicit symmetry breaking terms or coupling with ex- 
ternal fields. Initial investigation seem to reveal that the 



fractional states in the iV = 1 Landau level are direct 
analogues of the Laughlin state observed in GaAs. The 
1 /3 states in second Landau level are fairly robust unlike 
the case of GaAs, this is to be expected as the relativis- 
tic form factors for A = 1 lead to a more stable state 
when compared with the lowest Landau level 1/3-state 
in graphene.'^ 

As the sample mobilities increase it is expected 
that additional plateaus associated with the hierarchical 
states will also appear. A more interesting aspect would 
be the appearance of some non-Abelian FQH states such 
as the Pffafian or Read-Rezayi states.!!^ Unlike semicon- 
ducting 2DEGs, modification of interaction strength in 
graphene due to substrate manipulation could allow for 
tunability and hence stability of several FQH states. It 
has been proposed that modification of the dielectric 
environment could lead to stability of the Pffafian and 
Laughlin states.^ The chiral band structure of massive 
graphene and bilayer graphene also allows for the pos- 
sibility of tunable quantum phase transitions between 
incompressible and compressible states along with sub- 
tle transitions between incompressible states character- 
ized by different topological orders.I^ Furthermore since 
a metallic gate is used to vary the doping concentration 
an image charge accompanies the electron leading to elec- 
tric dipole-dipole interactions between the carriers which 
can be tuned by varying the distance of the the gate from 
the graphene layer. This dipolar interaction will also have 
important consequences on the edge states in particular 
with regards to the edg e re construction associated with 
FQH states of graphene.^i^oi 'pj^g important issue of FQH 
edge reconstruction has b een p reviously di scus sed both 
in the semiconductor 2DEGti^and graphen^^^llcontexts. 



VII. CONCLUSIONS AND OUTLOOK 

By demonstrating that a Landau level is pinned to 
the carrier neutrality point, quantum Hall studies of 
graphene's two dimensional electron system have in a 
very direct way confirmed the massless Dirac character 
of its bands. Similar studies have already also confirmed 
the massive Dirac character of bilayer graphene bands 
and are making progress toward revealing the distinct 
bands of ABA and ABC trilayers.!^ Indeed the few-layer 
family of graphene 2DEGs has impressive variety which 
can be accessed by varying the number of layers and how 
they are stacked. Electronic properties at low energies 
can depend on inter-layer coupling details that are not 
necessarily predictable with sufficient accuracy using ah 
initio electronic structure theory. We expect that sample 
characterization based on Shubnikov-de Haas oscillations 
and quantum Hall effects will continue to prove valuable 
in accurately characterizing the low-energy band struc- 
tures. Although this practical aspect of quantum Hall 
measurements is still important, studies of regimes in 
which interactions play an essential role are potentially 
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even more valuable because they may reveal interesting 
new physics. 

Quantum Hall studies have been shedding light on in- 
teraction physics of these systems in three different ways. 
First of all, the pattern of quantum Hall effects as the 
weak magnetic field limit is approached sheds important 
light on the zero-field ground state. This is particularly 
true of B = states that have sublattice polarization and 
associated momentum space Berry curvatures that lead 
to quantized anomalous Hall effects with Hall conductivi- 
ties (given in /h units) v = ±1, and v = ±2 in single 
layer graphene, and 0, ±2, and ±4 in bilayer graphene. 
When any of these states occur, the gap is expected to 
persist to finite fields and to occur at Landau level filling 
factor equal to the total Hall conductance integer which 
characterizes the B — state. As the field strengthens 
the transport phenomenology of the gapped state is ex- 
pected to cross over to that of an ordinary Hall plateau. 
If the ground state of single layer graphene was a CDW 
state, a, V = quantum Hall effect would be expected at 
all field strengths. The experimental observation that the 
I' = graphene quantum Hall effect does not survive to 
B — ?> 0, seems to establish unequivocally that the B = 
ground state is not a CDW. In bilayer graphene, on the 
other hand, many experimental hints are emerging that 
quantum Hall effects may survive to B = 0, although this 
question is just now in the midst of being clarified exper- 
imentally. Persistent v = ±4 quantum Hall effects at 
low fields suggest stable quantized anomalous Hall states 
whereas persistent i' — states at low fields suggest anti- 
ferromagnetic states with opposite spin-polarizations on 
opposite layers. The quantum Hall physics of graphene 
trilayers, which we have not reviewed extensively here, 
is still at an earlier stage of exploration. It is neverthe- 
less clear that ordered electronic states are likely even at 
B = 0, especially in ABC trilayers, and that anomalously 
persistent weak-field quantum Hall effects will if present 
play a valuable role in sorting out their character. 

Interactions in C2DEGs also play a prominent role by 
generating additional integer quantum Hall plateaus that 
are not expected by Landau quantization alone but are 
unrelated to possible B = spontaneous Hall states. 
Instead they can be understood by considering interac- 
tions among electrons that share the same well-resolved 
Landau level. Recent experiments on graphene sam- 
ples on h-BN substrates seem to suggest that interac- 
tions within the four-fold degenerate Landau levels lead 
to the formation of interaction induced gaps and hence 
additional quantum Hall effects. If these states spon- 
taneously break the SU(4) symmetry then the possibil- 
ity of observing exotic topological excitation associated 
with SU(4) Skyrmions remains intriguing. Recent exper- 
imental studies indicate that the v — QH plateau in 
graphene is an insulating state, indicating that it is most 
likely some sort of a density wave ordered state. More un- 
usual states are likely to be associated with Hall plateaus 
that appear within the = Landau levels of C2DEGs 
with a chirality index J > 2. In these cases the zero- 



energy Landau level supports multiple orbital degrees 
of freedom along with the four-fold degeneracy already 
present due to spin/valley. Speaking loosely, quantum 
states corresponding to cyclotron orbits with different ra- 
dius, which would have different energies in an ordinary 
two-dimensional electron gas, are degenerate. Experi- 
ments have confirmed the role of interactions in bilayer 
graphene as additional integer quantum Hall plateaus 
due to interaction induced gaps within the eight-fold de- 
generate zero-energy Landau level have been observed 
by several groups. When the Landau level orbital degree 
of freedom is active (for example at odd integer filling 
factors in J = 2 bilayers) Landau level orbital ordering 
reveals exotic QH states, which in the pseudospin lan- 
guage are analogs of spiral and helical ferromagnets. Be- 
cause the quantum Hall effects at odd integer filling fac- 
tors are least strong and more easily altered by disorder, 
they have not yet been studied extensively. Transport ex- 
periments in separately contacted bilayer graphene, and 
optical experiments are likely to yield many surprises 
in the future. Initial theoretical investigations in ABC- 
stacked trilayer graphene, which at the relevant energy 
scale belongs to the class of J = 3 C2DEGs, also indicate 
that integer induced quantum Hall plateaus which sup- 
port one-dimensional charge density waves with broken 
translational symmetry appear at certain integer filling 
factors.Iifii These type of charge density wave states at 
integer filling factors have no direct analogue in semi- 
conducting 2DEGs and are a consequence of interactions 
between electrons in different Landau level orbitals. The 
study of quantum Hall ferromagnets in bilayers and tri- 
layers is still in its infancy and appears capable of reveal- 
ing new physics that has not yet been anticipated. 

Finally of course, interactions in graphene systems of 
sufficiently high quality are expected to give rise to quan- 
tum Hall effects at a large number of fractional filling 
factors. Recent improvements in graphene sample mobil- 
ity have enabled the first observations of FQH states in 
graphene. The FQH states observed so far in the A^ = 
Landau level are most likely multicomponent generaliza- 
tions of the Laughlin states that are prominent in ordi- 
nary two-dimensional electron gases. The competition 
between incompressible states and disorder is certainly 
influenced by the difference in Landau level degeneracy, 
and this could be important in explaining the fractions 
which have been observed to date. As expected, due to 
the orbital character of the graphene which influences the 
Haldane pseudopotentials, the FQH states in the A^ = 1 
Landau level are m ore r obust than the FQH states in 
A^ = Landau level.1221221 The A^ = 1 FQHE in graphene 
is different from A^ = 0, and also from n = 1 in an or- 
dinary 2DEG, because of differences in the interaction 
form factors. An especially exciting possibility is the ob- 
servation of non-Abelian FQH states in graphene, which 
could end up being more robust than they are in the ordi- 
nary 2DEG case due to the advantage of tunable interac- 
tio ns in C2DEGs, as suggested by recent theoretical stud- 
^gg | 98 | 99 | Both the filling factors at which incompressible 
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states appear, and the character of their wavefunctions 
can change as a direct result of the chiral nature of charge 
carriers in graphene and bilayer graphene. Searches for 
correlated FQH states in the N — Landau level of bi- 
layer graphene with its additional orbital degeneracy is 
quite possibly the most exciting avenue and very likely 
to yield many surprises. Progress to date has been ham- 



pered by the numerical challenges posed by the additional 
degeneracies, and by the lack of predictive theoretical in- 
sights. 
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